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1 INTRODUCTION 



In this paper, we consider the dynamical problem of two rigid bod- 
ies interacting mutually according to the Newtonian laws of gravi- 
tation. It is well known that if both these bodies are approximated 
by point-masses, their dynamics are analytically solvable in terms 
of the integrable Kepler problem. However, even if only one of 
the bodies is an arbitrary rigid body then the complexity of the 
model changes dramatically. First of all, we have three more de- 
grees of freedom required to describe the rotational motion of the 
rigid body. Moreover, because the gravitational potential of a rigid 
body, in general case, depends on an infinite number of parame- 
ters, the two rigid-body dynamics depend also on this infinite set of 
parameters. Finally, orbital and rotational degrees of freedom are 
coupled non-linearly, and the resulting system is not analytically 
solvable. 

Because of the mathematical complexity involved, many ap- 
proximate models were studied in the past. For example, under cer- 
tain assumptions, one can consider the co-planar two rigid-body 
problem, see e.g., |Belets kij| (T975); | Barkin| ([1985}; [Beletskij &] 
|Ponomareva| ( | 1990) 1; Gozdziewski ( 2003jl. Another way to simplify 
the problem relies on truncating series expansions of the gravita- 
tional potential of the rigid body (the potential of mutual interac- 
tions), see e.g., |Markeev] ( |1967b|a| |1985 1988|; |Shcherbina| ( [T989) . 
In the so-called satellite (or restricted) models, it is assumed that the 
rotational motion does not influence the relative, Keplerian orbit of 
the mass centres of the bodies. In such models, the relative orbit 
is given parametrically, and only the rotational dynamics need to 
be investigated; for details, see e.g., Markeev ( 1985 196 7bl[T975^ ; 
|Maciejewski (1994) ; |Maciejewski & Gozdziewski (1995 i. 

In this work, we focus on particular version of the model in- 
volving a sphere (a point-mass), by considering that the rigid body 
is axially symmetric. Obviously, this leads to much simplified form 
of the gravitational potential. However, in general, it still depends 
on an infinite number of parameters. For the first time this problem 
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was investigated by |Kinoshita| ( |1970 ^, and hence it will be called 
the Kinoshita problem from hereafter. Our aim is to determine all 
possible classes of the relative equilibria (stationary solutions) in 
this problem, using a general formalism in [Wang et al.|p99T| and 
|Maciejewski| \ 1 995\ . Keeping in mind a motivation of the work by 
Kinoshita, we try to generalize his analysis and to simplify condi- 
tions for the existence of these equilibria, avoid repetition of "of 
assuming only the necessary general form of the gravitational po- 
tential of the rigid body. In this sense, our paper is also related to the 
work by |ScheereS| ( [2006^ , who searched for the relative equilibria in 
the unrestricted two rigid-body problem, and formulated conditions 
of the equilibria in terms of solutions to two non-linear equations, 
parameterized by integrals of the total energy and angular momen- 
tum. Moreover, we stress that our analysis rely on the very basic 
vectorial form of the equations of motion, and we consider the full, 
non-restricted model with symmetric rigid body. The assumption of 
symmetry implies that the results derived for the general, full two- 
body model cannot be simply "translated" as a particular case, and, 
in fact, a special reduction of the system which takes into account 
the symmetry explicitly is necessary. We note that the relative equi- 
libria of the model with axially symmetric body may be understood 
as particular periodic orbits in the full two rigid body model con- 
sidered by Scheeres^ (200 6 1. This reduction is in fact crucial for a 
generalization of the results of Kinoshita and to discover a class of 
equilibria that have been missed in his analysis. 



Investigations of the unrestricted two rigid body problem, in 
a version considered in this work, and its particular stationary so- 
lutions, are interesting because they concern a special case of gen- 
erally unsolved, classic problem of the dynamics. Moreover, the 
qualitative analysis of this model may help us to answer for impor- 
tant, even "practical" questions on the dynamics and coupling of 
the rotational and orbital motions in many astronomical systems. 
These are, for instance, binary asteroids, see |Scheeres| j2006] > and 
references therein. A deep understanding of the qualitative dynam- 
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ics are important for the attitude determination and control of large 
artificial satellites orbiting planets and/or irregular natural moons. 
Recently, the diversity of extrasolar planets discovered in wide dy- 
namical environments also rises questions on their rotational mo- 
tion and related long term effects (e.g., |Correia et al.|2008[ >. In the 
context of mathematical complexity of the problem, the relative 
equilibria are the simplest class of solutions that may be found and 
investigated analytically, and are helpful to construct local, precise 
analytic theories of motion in their vicinity in the phase space (e.g., 
[Gozdziewski & Maciejevyski|1998| ). 

A plan of this paper is the following. In Sect. 2 we formulate 
the mathematical model and the equations of motion in the most 
general, vectorial form are derived. Section 3 is devoted to define 
the Kinoshita problem. Next, we introduce the relative equilibria, 
and we perform a global analysis of their existence and bifurcations 
(Sect. 4). This part relies on particular, geometric reduction of the 
equations of motion, and is the primary key for our analysis. Fi- 
nally, we compare the classic results by Kinoshita with the results 
obtained through the approach introduced in this work. 



2 THE EQUATIONS OF MOTION 

Let us consider the gravitational two rigid body problem. We as- 
sume that one of the bodies, CP is a point mass, or an uniform sphere, 
with mass mi. The second one, 25, is a rigid body with mass m2 
(Fig.[TJ. The mechanical problem in the most general settings has 
nine degrees of freedom, so the dynamics of the bodies are de- 
scribed by 18 first-order differential equations. Moreover, the dy- 
namical system possesses symmetries, i.e., it may be shown that the 
equations of motion are invariant with respect to a six-dimensional 
group. This fact is a direct consequence of the very basic laws of 
Newtonian mechanics. Namely, the equations of motion do not de- 
pend on particular choice of the inertial reference frame. Thus, we 
can choose the origin of the inertial frame in an arbitrary point, and 
this implies that the equations of motion are invariant with respect 
to the natural action of three dimensional group of translations. 
Moreover, also the orientation of the inertial frame can be cho- 
sen arbitrarily. Hence, the equations of motion are invariant with 
respect to the natural action of the three-dimensional group of ro- 
tations. 

The existence of these symmetries makes it possible to reduce 
the dimension of the phase space of this system and to simplify 
its analysis. However, there is no unique procedure for such a re- 
duction. Obviously, it should rely on such a transformation of the 
initial system of the equations of motion (phase-space variables) 
that the resulting equations form a lower-dimensional sub-system. 
In fact, the reduced equations of motion of the two rigid body prob- 
lem were obtained by Wang et al.| ( [T99T| l (actually, they considered 
the motion of a rigid body in the central gravitational field), and, 
with a simpler and more direct method and under general settings 
by |Gozdziewski & Maci^ewski| ( |1999) . These papers are good ref- 
erences to step-by-step development of the reduced equations of 
motion, which is skipped here to save space. 

Before we write down these equations of motion, we need 
to fix the notation. Components of a vector x (as a geometric 
object) in an inertial reference frame will be denoted by x = 
[ji,X2,X3]^. Components of the same vector in the rigid body 
fixed frame we will denote by the corresponding capital letter, i.e., 
X = [X\ ,X2,X2]^ . Thus, if A is the orientation (attitude) matrix of 
the body with respect to the inertial frame, then we have x = AX. 
The scalar product of two vectors x and y is denoted hy x-y. It can 




Figure 1. Geometry of the problem. X, Y, Z are units vector coincided with 
the principal axis of inertia of the rigid body. 



be calculated in an arbitrary orthonormal frame as follows: 

3 3 
x y = (x,y) := Y^x^ji = (AX,AY) = (X, Y) = £x,-y;. 



We shall also write: 



--x-x^ {x,x) =X^ ■- (X,X). 



In our exposition, we follow [Maciejewskil j 1 995^ ; |Gozdziewski &| 
|Maciejewski| ( [1999l l. The reduced equations of motion describe the 
relative motion of the bodies with respect to a frame fixed in rigid 
body S. They have the following form of the so-called Newton- 
Euler equations: 

d 

— R = Rxf2 + P, 



At 3R 



At 



G = Gxf2 + Rx 



dU_ 

5r 



(1) 



where R is the radius vector directed from the mass centre of S to 
7, P is the relative linear momentum, G is the angular momentum 
of 2, and fl = I^'G is its angular velocity; U is the gravitational 
potential of the body 23, and I = diag(A, S,C) stands for its tensor 
of inertia. We assume that the rigid body fixed frame coincides with 
the principal axes frame. The geometry of the model is illustrated 

in Fig. [T] 

As it was shown in |Gozdziewski & Maciejewski| ( |1999} , equa- 
tions l|T]( are Hamiltonian with respect to a non-canonical Poisson 
bracket. They admit two first integrals: the energy integral. 



«=^(P:P) + ^(G,£^)+(/(R), 



(2) 



and the modulus of the total angular momentum L, which, with 
respect to the body frame, is given by: 



l2 = (L,L), where L:=RxP + G. 



(3) 
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The total angular momentum with respect to the inertial frame, i.e., 

I = AL, (4) 

is a constant vector. The time evolution of the attitude matrix A is 
given by the following kinematic equations: 



(5) 



where, as i n|Gozdziewski & Maciejewski| ( |1999] l, for a vector X - 
\X\ , X2 , X^Y , we denote: 



X:= 






-X3 


X2 


^3 









^1 






(6) 



If the rigid body 23 does not possesses any additional symmetry, 
then the relative equilibria of the considered two rigid body prob- 
lem are just the equilibria solutions (critical points) of the reduced 
equations l[TJ. They were investigated by many authors, see |Wang] 
[eTa l. ( 1991 ); Scheeres (20061 and references therein. The most in- 
triguing result of these investigations is the discovery of the so- 
called non-great-circle, or non-Lagrangian solutions. If the system 
is in a relative equilibrium, then the relative orbit is circular. How- 
ever, in a case of non-Lagrangian equilibrium, the orbital planes 
of the point ? and the rigid body S are mutually parallel but do 
not coincide. These classes of equilibria have been discovered by 
Abulnaga & Barkm] jl979[ l. Later on, they were re-discovered by 



WangetaL| ( |i99T] r 



3 THE SYMMETRIC KINOSHITA PROBLEM 

From hereafter, we shall assume that the rigid body B is dynami- 
cally symmetric with respect to its third principal axis, i.e., we set 
I = diag(A,A,C). Moreover, we assume also that the gravity field 
of 2 is symmetric with respect to the same axis and that it can be 
expressed through the following condition: 

dU 



(R) 



:0. 



(7) 



(We note, for a further reference, that the equatorial plane of the 
rigid body is defined in such a way, that it contains the mass centre 
and this body and the symmetry axis is normal to that plane). Under 
these assumptions, equations (|TJ have an additional first integral 
G3. The existence of this first integral is connected with the fact 
that now the reduced equations l[TJ have one-dimensional symme- 
try group. To show this explicitly, we need to introduce additional 
notations and to state or recall some elementary facts. 

By SO(3,R) we denote the special orthogonal group, i.e., the 
group of 3 X 3 matrices A, satisfying AA^ = E, and detA = 1. Its 
Lie algebra so(3,K) consists of 3 x 3 antisymmetric matrices. It 
can be identified with R^ by the following isomorphism: 



3Xh 



■X:= 



-Xi X2 
X3 -Xi 
-X2 Xi 



e so(3, 



the 



Later on, we use several identities which we collect as the following 
statements (e.g., |Maciejewski|1995[ l. 

Proposition 3.1. For arbitrary X, Y e R^ and A e S0(3 
following identities are hold true: 

(i) X^ = -X = {^). 

(ii) XX = XxY = -YX. 



(iii) XY^ YX^-X^YE. 

(iv) XxY = XY - YX. 

(v) AX = AXA^. 

(vi) A(X x Y) = AXY = (AX) x (AY), 
where E denotes the 3x3 identity matrix. 

All rotation about the third axis of the body form a subgroup of 
SO(3,R) isomorphic with SO(2,R). We show that equations (T) 
possess SO(2,R) symmetry. That is, they do not change its form 
after the following transformation: 

[R,P,G] = [A3R*,A3P*,A3G*1, (8) 

where A3 is an arbitrary matrix of rotation about the third axis of 
the body, i.e., a matrix of the following form: 

-b 



A3: 



2 2 
where a -\-b = I. 



(9) 



Let us notice that A3I = IA3 (i.e.,I commutes with A3), hence we 
have, 

^R = A3^R* = (A3R*) X (r'A3G*)+A3P*. 



xI-'G* 



This equation may be rewritten as 

(A3R*) X (r' A3G*) = (A3R*) X (A3r 'G*) = A3 [r* 

where the last equality follows from the sixth identity given in 
Proposition |3.1| Hence, 

A3 — R* = A3 [R* X fj* + P*] , 



where Q.* 



I 



'dt 

^ ' G* , and, finally, 
d 



-R* =R* xfl* 



dt 

and this shows that the first equation in l[T]l is invariant with respect 
to variables change ([Sj. In a similar manner, we show that the re- 
maining two equations have the same property. 



4 RELATIVE EQUILIBRIA 

4.1 General tlieory of the relative equilibria 

Let us recall a formal, geometric definition of a relative equilib- 
rium of a system with symmetry (see, for instance [Marsden &| 
|Ratiu|[l"994) , which is crucial to perform our analysis: a relative 
equilibrium is a solution of the system represented by the phase- 
space curve which is an orbit of a point under the action of a one- 
dimensional subgroup of the symmetry group of the system. 

In our case, the symmetry group is S0(2,R) which is identi- 
fied with matrices of the form It is a one-dimensional group. 
Let us put a = cos(p and = sincp in equation l|9j. With this param- 
eterization, we denote elements of S0(2,R) by A3((p). Hence, an 
orbit of a point [Ro,Po,Go]^ under the action of SO(2,R) is the 
following set: 

{[A3((p)Ro,A3((p)Po,A3(9)Go] gR"* | (pe [0,27r)}, 

where Ro,Po,Go are constant vectors. Thus, the relative equilib- 
rium of the symmetric Kinoshita problem is a solution to equa- 
tions l|T](, which has the following general form: 

R(0=A3(M)Ro, 

P(0=A3(M)Po, ) (10) 
G(0=A3(M)Go, 
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where N is the real number Notice, that £l{t) = I^'G(/) = 
A3(A';)f2o' where flo = I^'Gq. In other words, the relative equi- 
librium is a special periodic solution to the equations of motion l[TJ. 
From the above formulae we can deduce more useful conclusions 
on the relative equilibrium: 

(i) vectors R(/), P(f), G{t) and Q.{t) have constant lengths, 

(ii) their "third" components R^it), /'3(f), G^it) and Q.2{t) are 
constant, 

(iii) all angles between vectors R(?)> P(f). G(f) and Q.{t) are 
constant. 

Let us notice that the relative equilibrium \10\ does not determine 
an unique set of [Rq , Po , Go] . In fact, instead of [Rq , Fq , Go] , we can 
take [A((p)Ro, A((p)Po, A((p)Go], with arbitrary chosen constant (p. 
One would like to perform a reduction of l|T]l in such a way that 
in the reduced system, the relative equilibrium corresponds to an 
"usual" equilibrium point (i.e., which is understood as the criti- 
cal point of the equations of motion). It is possible, however, it 
may lead to singularities. For example, let us introduce three sets 
of cylindrical coordinates: 



R= (pcos(p,psin(p,/?3)"^ , 
P= (pcos-u,psin\),/'3)^ , 
O. = (cocos\|/,(Bsin\|/,tl3)^ . 



Using these coordinates, we may write: 



^ = (L/pCos(p,f/pCos(p,(/3) , 



where. 



t/p = ^, and f/3 = ^. 

With respect to these new variables, system (|T]l reads as follows: 
d 



df 
d 
At 
d 
At 
d 
At 

d 

df 
d 
At 
A 
At 
A 



p = pcos('u — (p) — zcosin(\|/ — (p), 
z = P3+pa)sin(\|f-(p), 

v> = [/'3a)cos(\(/— (p) + (/psin(\|/— (p) — p£lo, 



-/'3C0sin(\|/ — (p) — Up cos(t) — (p), 
Pt, = pcosin(\|/ — u) — U-i, 
00 = A"' [zUp-pUi] sin(\|/-(p), 
9 = P ' [/-?sin(-u - (p) +za)cos(\|/- (p) - pO.^] , 



-\|/=(Aco)-i [(z(/p-pC/3)cos(\|/-(p)-(A-C)con3] , 



Hence, finally we have one equation less than in ([TJ because Q.j, is 
constant. Note that the right hand sides of the equations depend on 
the difference between the angles but not on the angles themselves. 
Thus we can introduce the following new set of variables: 



and then we obtain: 
d 
At 
A 
At 
A 



-z = 



At 



At 
A 

At' 



At 



/7Cos((3) — zcosin(a), 

P3 + pa)sin(a), 

p = -P3a)sin(a-p) - f/pCos(P), 

pcosin(a-|3)-(/3, 

-co = \zUp - pUi] sin(a), 

zUp - PU3) CO"' cos (a) + (Cf23- 
Aco) [psin(p) +zcocos(a)] p"', 
P = p"' [/'3C0cos(a-P) + (/psin(P)] - 
p"' [psin(P)+zcocos(a)], 



d 

d7'P = «- 

The first seven equations in the above set represent a system of 
equations which flows from the reduction l[T]l by the symmetry. 
Now, a relative equilibrium of ([TJ corresponds to an equilibrium 
point of the above system. However, this system has singularities 
when one of variables p, CO or p vanishes. Because of this fact it 
would be very difficult to perform global, qualitative analysis which 
is our primary goal. 

4.2 The geometric reduction of the system 

A possibility that the reduction may introduce singularities to the 
equations of motion inspired us and, in fact, forced us to choose 
a particular approach which we describe below in detail. Clearly, 
instead of using any local variables, it is much more convenient to 
derive the necessary and sufficient conditions for the existence of 
relative equilibria in terms of the original global and non-singular 
variables (R,P,G). 

Let us assume that a relative equilibrium is given by |TOj. Our 
aim is to find equations determining (Ro,Po,Go), and A'. At first, 
we may notice that: 



— A3(M) 
At ' 



A3(M)N, 

where N = [0,0, A']^. Thus we may easily derive that: 

A3(M)NRo = A3(M)(N x Ro) 



'^-RM = ^A3(A'f)Ro^ 



At At 
On the other hand, we have also: 
d 



At 



R(r) = R(f) X a{t) + F{t) = Ai{Nt)(Ro x fio +Po) 



as the system has S0(2,R) symmetry. So, we have 
NxRo = Roxno + Po. 

We proceed in a similar way with other components of the 
phase variables defining the relative equilibrium, i.e., with P(;) and 
G(t). Finally, we obtain the following equations: 



N X Ro : 
N X Po : 



-- Rq X flo + 

: Po X flo - 



Po, 

aRo' 



a = \|/ — (p, 



, = 'U-(p, 



Nx Go = Go X flo + Ro X 
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This system can be rewritten in a more compact form as follows 

fixRo = Po, 

dU 

aRo' 



n X Po : 



:Rnx 



9Ro 



(11) 



where, H = fio + N, is the instantaneous angular velocity of the 
rigid body. 

System \\\\ defines nine equations which must be satisfied by 
ten unknowns. We found that there are several ways to overcome 
this inconvenience. In our further procedure, we choose the follow- 
ing approach. We look for the relative equilibria which exist for a 
given value of the angular momentum of the system, chosen as the 
model parameter. Because 

l2 = (R X P + G)2 = (R() X Po + Go)^ , 

after the elimination of Pq, we obtain: 

ai?g-Ro(n,Ro)+Go 

Hence, our goal is to find all solutions to the following nonlinear 
equations: 

OxRo =Po, 

dU 



i2(i2,Ro) -Rofi = - 



aRo' 

du 

i2xGo = R()X 3—, 

dKo 



(12a) 
(12b) 



fli?,^-Ro(fl,Ro)+Go) =l}. 



(12c) 



(12d) 



In the above equations, L, and the third component of Go are now 
fixed parameters of the reduced model. 

4.3 Geometric description of the equilibria solutions 

The equations of motion l[TJ describe the relative motion of the bod- 
ies in the rigid body fixed frame. Apparently, this makes it difficult 
to obtain the geometrical and physical interpretation. In fact, we 
can deduce all necessary information quite easily. 

First of all, let us recall that the total angular momentum / is a 
constant vector. So, we can choose an inertial frame in such a way 
that its third axis is directed along that vector. If R(f) = A3(A'?)Ro 
describes the relative orbit in a relative equilibrium, and r{t) is the 
corresponding equilibrium vector, then we have: 

l-r(t) = (l,r(/)) = (AL,AA3(M)Ro) = (L, AalA'ORo). 

Moreover, the following general relation holds true: 

L = A3(iVO(RoxPo + Go), 

so, finally, we have: 

l-r{t) = (Ro xPo + Go,Ro) = (Go,Ro). (13) 

We can conclude that in a relative equilibrium, the projection of 
the relative radius vector onto the total angular momentum vector 
is constant. An important conclusion follows immediately: if this 
projection is non-zero, then orbits of the point and the rigid body 
lie in different planes. 

Letz(f) denotes the unit vector along the symmetry axis of the 
body. Then 

l-z{t) = (L,Z) = (A3(M)(Ro X Po + Go),Z) = (R,, x Po + Go,Z), 



because Z = [0,0, 1]^, and so, Ki{Nt)Z = Z. Thus, in a relative 
equilibrium, the projection of the symmetry axis onto the total an- 
gular momentum is also constant. 

Let us introduce the orbital reference frame with axes parallel 
to r{t), r{t) and c{t) := r{t) x r(t). Because the relative orbit in a 
relative equilibrium is circular, this frame is orthogonal. We show 
that the projection of the symmetry axis z{t) onto the axis of this 
frame is also constant. In other words, we are going to prove that 
in the relative equilibrium the symmetry axis of the rigid body has 
a fixed orientation with respect to the orbital reference frame. 

To shorten notation, we shall write A3 instead of A^iNt). Let 
B = AA3. Then we have: 



df 



B = AA3 + AA3 = Af2A3 + AA3N. 



But in a relative equilibrium we have = A3f2o- By the fifth prop- 
erty given in Proposition |3.1[ 



hence we obtain that: 



B = Bi2, 



df 

where fl = flo +N is a constant vector. Now, it is easy to show that: 
r(0=AR(?)=BRo, 

r(0=B(iixRo), (14) 
r(Oxr(0=B(Rox(fixRo)). 
Because z{t) — AZ — BZ, and A3Z = Z, we have: 

r{t)-z{t) = (r(0,z(f)) = (BRo,BZ) = (Ro,Z), 
r{t)-z(t)^{0.xRo,Z) = {Po,Z), (15) 
c{t) ■ z{t) = {Rox {a xRo),Z). 
The above formulae prove our initial claim. 

5 CONDITIONS FOR THE RELATIVE EQUILIBRIA 

In this section, we show that all relative equilibria can be deter- 
mined by solutions of certain non-linear system comprising of only 
two non-linear scalar equations. We shall also prove that there ex- 
ist three classes of these relative equilibria in the dynamical model 
which we consider. 

In order to simplify notation while working with coordinates, 
we shall skip index "0" when denoting coordinates of vectors Rq, 
Po, Go from hereafter. Thus, we set 

Ro = [Ri,R2,R3f , 

flo = [i2l,i22,i^3]^, 
0.= [aiM2,0.3+Nf , 
Go = [AQ.i,AQ.2,CQ.3f , 

5^:=3^(Ro) = [f/i,f/2.f/3f. 

The fact that the body is symmetric implies that: 

RiU2-R2Ui=0, (16) 
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with accord to |7}. Taking the vector product of both sides of equa- 
tion ( |12b| ( with R(), we obtain: 



(i2,Ro)iJxRo = Rox 



dU 

aRo' 



(17) 



Hence, talcing into account l |16^ , we have: 

(n,Ro)(Z,OxRo)=0. 

Thus, we have to consider two different cases. 

If (f2,Ro) = then from l |17| l it follows immediately that: 

or, (Z.fi X Ro) = 0, and this gives 

Q.iR2-^2Ri =0. 

In the first case, the gradient of U at point Rq is parallel to the radius 
vector. According to |Scheeres| ( |2006| l, such a point is called locally 
central. We shall say that a relative equilibrium is locally central if 
point Ro is locally central. 

In our analysis, we assume that the gradient of the potential 
does not vanish at a considered point Rq. Thus, a relative equilib- 
rium is locally central if and only if (O.Rq) = 0. In fact, if Rq is 
locally central, then from it follows that either {Q., Rq) = 0, or 
n X Rq = 0. If the second possibility occurs, then, by (Ha) , Pq = 0, 
and, as equations \\\) imply, the gradient of U at Rq vanishes. But, 
according with our assumptions, it is impossible. 

At the final stage of our analysis, it is convenient to use cylin- 
drical coordinates (p, 'SfyRj,) instead of [Ri ,R2,R-i). The potential U 
expressed with respect to the cylindrical coordinates depends only 
on p, and Ry, still, as a function of these two arguments it will be 
denoted by the same symbol U . It will be clear from the context 
which coordinates are in use. 



5.1 Locally central relative equilibria 

If Ro is locally central, then, as we have shown, (O.Rq) = 0, and 
of course: 

«''^aRo=»- 

Thanks to this fact, the basic system of the equations of motion \\2) 
may be simplified considerably, and it reads as follows: 



Po = nxRo, 
X Go = 0, 



ni?^ + Go) =L?-. 



(18a) 
(18b) 
(18c) 
(18d) 



Equation (|18c^ is equivalent to the following two scalar equations: 



\AQ.j,-CO.^ \ 0.2 = 0, 



\AQ.j,-CQ.^ \ Q.\ =0. 



5.1.1 Case Q.\ = Q.2 = (cylindrical precession) 

\f Sl\ =£l2=Q then condition (fl.Ro) = reduces to £l^R^ = 0. 
But £^3 7^ 0, otherwise O.^ = 0, and, as l |18b[ l shows, the gradient 
of U vanishes and then Rj = 0. Moreover, Ri and R2 do not van- 
ish simultaneously, otherwise Ro = 0. Hence, we can safely use 
cylindrical coordinates (p,z = o|^because p > 0. We introduce two 
functions of p given by: 

/(p):=(/3(p,0), and g{p) := Up{p,0). 

Now, equations (fTSbll and (TTSdl lead to the following system: 



/(p) = 0, pm=8{p) 



[p^n3+c£^3 



(19) 



Through appropriate simplification, we find that the above system 
is equivalent to the following ones: 

/(p) = o, {L+Eca,f = p'g{p), a, = -e^^±^, 

(20) 

where parameter e e {—1,1}. 

The above considerations can be summarized in the following 
way. In the considered case, the relative equilibrium exists if and 
only if there exists a solution p > of equations: 



/(p) = 0, (L + em3)' = p3^(p). 
If such a solution exists then we define: 



i?i=pcos(p, /?2 = psin(p, a^ = p-^(/p, 



(21) 



where (p e [0, 2n] can be chosen arbitrarily. Then the relative equi- 
librium is determined through: 



Rn 



"^1" 




-R2 




■ " 


Ri 


, Po = P 


Ri 


, Go = C 















Q.3 



From the above formulae it immediately follows that point 7 and 
the mass centre of the rigid body 25 move in one plane. In fact, 
by IT3), we have: 

/•r(r) = (Go,Ro)=0. 

Moreover, by ijTS), we have also: 

r{t)-z{t) = {Ro,Z}=0. 
i-(0-z(f) = (Po,Z)=O. 

Hence, the axis of symmetry of the rigid body is perpendicular to 
the plane of the relative orbit, see Fig.|2] We called this kind of the 
relative equilibrium as the cylindrical precession from hereafter. In 
the paper of Kinoshita, it is named, after |Duboshin| l 1959), as the 
float solution. Regarding existence conditions of this solution, let us 
notice that system {21) is overdetermined, we have two equations 
for one variable p. However, if the body is symmetric with respect 
to the equatorial plane, then U is an even function of i?3. In this 
instance, we have: 



m- 



du 



(p,0)=0, 



identically. Thus, in this case, the cylindrical precession exists if 
and only if the second equation of system \2l\ has a solution. 
Moreover, if the potential is such that there exists p > 0, for which 



Therefore, either Sli = Q.2 = 0, or Afl3 — CQ.T, = 0. We investigate 
both cases separately. 



' Note that due to the symmetry, the third angular coordinate (|) is irrelevant 
here. 
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Figure 2. Cylindrical precession in the inertial reference frame. Shaded 
bodies show the attitude of the rigid body 23 and the point CP in other mo- 
ments of time. 

^(p) > 0, then we always may find parameters of the problem such 
that equations \2\) are fulfilled. 

If the body is not symmetric with respect to the equatorial 
plane, the first equation of system l |21^ is not satisfied identically. 
Thus a question emerges: does there exist an axially symmetric 
body which is not symmetric with respect to the equatorial plane, 
and for which equation: 

(/3(p,0)=0, 

has a solution p > 0? To answer this question, let us consider an 
axisymmetric rigid body composed of three uniform spheres whose 
centers of masses are placed at the same line. The masses of these 
spheres are chosen in a way to guarantee that the potential is not 
symmetric with respect to the equatorial plane. The resulting func- 
tion /(p) is depicted in Fig. |3] Since its graph crosses the p-axis 
for p > 0, there exist such points which satisfy the first condition of 
system ip-i) , indeed. Then the parameters of the system can always 
be chosen to fulfill the last condition of \2\) . Thus, the cylindrical 
precessions may exist in the system with the rigid body, which is 
not symmetric with respect to its equatorial plane. 

5.1.2 Case AQ.^ = CQ.^ ( inclined co-planar precession) 
Condition AQ.J, = CQ.T, implies immediately that 

A 

Moreover, now (ft, Rq) = 0, is equivalent to (Go, Rq) = 0, because 
Af2 = Go- We have even more, because now equation l |18b| ) can be 
rewritten in the following form: 

Po = ^GoxRo. (22) 
and equation l |18d| l reads as: 

l2 = (aRl + Gof = ^{rI+A)^. (23) 
Therefore, equation ijTSb} can be rewritten as follows: 

A solution of this equation gives Rq. However, this form is not con- 
venient for further analysis because solutions to the above equation. 
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Figure 3. The plot of /(p) crosses p-axis for p ~ 0.26. The values of 
parameters are chosen in the following way: m = 0.15, = 0.3. ,U3 = 0.55, 
where jjj = mi/Y,^j^i ™; rnass of the i'-th sphere. See the text for 

details. 

whenever they exist, are not isolated. In terms of the cylindrical co- 
ordinates, we may rewrite ( |24[ l as follows: 

P — ; 9 = Up, 

(Rj + p^+Af 

(W2 + p2+A)' 

A solution to the above system gives us Rq. Because, due to the 
symmetry, we may choose freely the polar angle, then we can also 
assume, without any loss of generality that R2 = 0. Vector Go is 
perpendicular to Rq, and its modulus is fixed by l |23| l. In the generic 
case, R\R3 7^ 0, and we can put: 

Go = ai A [0, 1 , 0] ^ + fl2 ^ [/^3 , 0, 1 ] ^ , 

where aj and 02 ars arbitrary real numbers satisfying the following 
relation: 

2 2 

(71 +(29 = T -. 

' ' [Rl+Ar 

For each choice of ai and a2, we have Go, and then the correspond- 
ing value of Pq is given by l |22^ . Thus, a solution to equations l |25| l 
gives us a one-parameter family of the relative equilibria. 

From the above formulae it immediately follows that point "P, 
and the mass centre of the body S move in one plane. In fact, we 
have: 

/•r(f) = (Go,Ro)=0. 

Moreover, as in the first case of the cylindrical precession, we can 
determine the orientation of the symmetry axis using {15}. We have: 

r(/)-z(/) = (Ro,Z) =^^3, 
r{t)-z{t) = {Fo,Z)=P3 = -aiRu 
c{t)-z{t) = (Ro X (n X Ro),Z) = -a2RoRi. 
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Figure 4. The inclined co-planar precession in the inertial reference frame. 
Shaded meshes are for an illustration of the attitude and relative positions 
of the rigid body 3 and the point CP at different moments during the motion 
of the system. 



Thus, in general, the symmetry axis of the rigid body is inclined 
with respect to any axis of the orbital reference frame. We shall call 
this class of relative equilibria as the inclined, co-planar regular 
precession, and its geometry is illustrated in Fig.|4] 

5.2 Non-locally central case (conic precession) 

As it has been already established, the non-locally central equilibria 
are characterized by P3 = R20.\ —R\£i.2 = 0. It appears that in such 
an instance, instead of using system \\2\ , it is better to analyse 
equations in which vector P is not eliminated. Hence, we consider 
the following system of equations: 



X Ro = Po, 



nxP, 



■- 



du 



f2 X Go = Ro X 



du 

3Ro' 



(RoxPo + Go)=l2 



(26a) 
(26b) 

(26c) 
(26d) 



Our aim is to find all solutions of this system which are not locally 
central. Thus, we look for solutions for which none of variables p, 
i?3, ^3 vanishes. Hence, we can express Ui and U2 in terms of Up, 
i.e.. 



Ui 



Ri 



for ;=1,2. 



Since vanishes, from the first two components of equations \26b\ 
it follows that 



Pi 



Up Up 
-Rl^, and P2=R\ 



(27) 



fi3P f23P 

Knowing Pi and P2, from the first two components of equation 
\26&) , we find: 



Sli=Ri- 



for (=1,2. 



(28) 



pRiQ.2 

Substituting {27\ , (|28j into the third component of equation |26bJ, 
we obtain: 



{pUp+R3U3)0.l = U^ 



(29) 



Now, let us consider equation l |26c[ l. It is easy to see that the third 
component of this equation is fulfilled identically. Using (|28j, the 
first two components can be rewritten in the following form: 

Ri (a0.3-CQ.3^ UpUi =Ri {pUp +RiU3) {pU^-RiUp) h^, for ;= 1,2. 

As Ri and R2 do not vanish simultaneously, the above system is 
equivalent to just one equation: 

AQ.2-C0.3^ UpUi = {pUp +R3U3) {pU3-R3Up) Q.3. (30) 

Next, substituting expressions {27\ and l |28[ >, equation l[26dj is 
rewritten as: 

l2 .2 



^2^ [{A+Rj)U3 + pR3UpY ^ fpUp ^ 



PUP+R3U3 



Q.3 



(31) 



Finally, we use l |3Q[ ) to eliminate Q3 from equations ( |29[ l and \31) . 
Thus, eventually, we obtain the final subsystem that makes it pos- 
sible to determine the phase-space variables of the relative equilib- 
rium, in quite a compact form: 



(32) 



L'-UiK = [U^ + Ufj [AU3 +R3Kr , 
C^a^f/|A:= [AUpU3+K{R3Up-pU3)f , 



where we introduced a new variable: 

K = pUp+R3U3. 

A solution to the above conditions gives us Ro. The remaining equi- 
librium variables can be determined by equations ([27j, l|28}, and 
\29\ , respectively. 

For this solution we can find, after tedious simplifications, that 



/•r(0 = (Go,R(,)^ 



Up{AU3+R3K) 
U3K 



R3UP-PU3) 



According with our assumptions, p 7^ 0, and hence Up 7^ 0. Thus, 
the right hand side of the above equation vanishes if either R3Up — 
PU3 = 0, or AU3 + R3K = 0. In the first case, Rq is locally central 
point which we excluded from our considerations. In the second 
instance, from equation ([32j, we obtain that L = 0, so it is also 
excluded from our considerations. We conclude that non-locally 
central relative equilibrium is non-Lagrangian ones. 

In order to determine the orientation of the symmetry axis in 
the orbital frame which is specific for this type of relative equilibria, 
we can easily find that: 

r(0-zW = (Ro,Z)=i?3, 
r{t)-z{t)=0, 



c(0-zW = (Rox(i2xRo),Z) = 



0.3 



Thus, in the orbital reference frame, the axis of symmetry always 
lies in the plane formed by the relative position vector and the nor- 
mal to the orbital plane. We call this kind of the regular preces- 
sion as the conic precession from hereafter. Notice, that in the con- 
sidered case, vector c{t) is not perpendicular to the orbital planes. 
Moreover, we have: 



l-z{t) = (Ro X X Ro) + Go,Z) = 



pUp^ 
Q.3 



C£l3. 



Thus, the axis of symmetry is inclined to the orbital plane. The ge- 
ometry of this solution is presented in Fig. |5] Unlike the locally 
central cases, even for bodies with simple potentials, it is difficult 
to decide whether equations {i2\ have a solution. In order to jus- 
tify that the relative equilibria belonging to the conic precession 
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Figure 5. Conic precession in the iiiertial reference frame. Shaded bodies 
show the attitude of the rigid body 2 and the point CP in other time points. 



does not vanish identically. Hence, by the implicit function theo- 
rem, system \32\ has a solution p — p(Q.2), R3 = R^iO.^), such that 
p(0) = 0, and R^iO) = Ry Moreover, such solution is unique. 

We underline an important fact that, even if the body is sym- 
metric with respect to the equatorial plane, the described conic pre- 
cession does exist. In other words, splitting of the orbital planes 
can be induced not only by an asymmetric mass distribution, it can 
be also induced by a proper, particular rotation of the rigid body. 
Remarkably, the first type of solutions has been known since |Ab^ 
|ulnaga & BarkinH19 79'>, who constructed them for a class of non- 
symmetric potentials and further analysed their consequences in the 
motion model of the Earth-Moon system l |Barkin|1980^ . 



6 A COMPARISON WITH THE WORK OF KINOSHITA 

Our analysis presented in the previous sections lead to somewhat 
different conclusions than those ones in ( |Kinoshita|1970l >, and m 
this section we give an explanation and overview of these discrep- 
ancies. Kinoshita found three classes of stationary solutions: 



class may really exist, we apply the following reasoning. If p, as 
well as ^3, tends to zero, then the conic precession tends to a par- 
ticular case of the inclined planar regular precession described in 
Section 5.1.2 In fact, if p = f23 = 0, then Up = 0, so the second 
equation in \32\ is satisfied identically, and the first one coincides 
exactly with the equation describing corresponding case of the in- 
clined co-planar solution. Thus, we can take this particular solu- 
tion as the zero-th order approximation (with respect to Q.^) of the 
conic precession solution, and we apply the perturbative approach 
to solve the system j32[ >. 

Thus, the zero-th order solution is given through: 

p = 0, R3=Ri, 
where R^ is a solution of the following equation 



U3 



p=o- 



[A+R: 



To shorten the notation, we define: 



U3 



Taking into account that 

f^p|p=()= t^p3|p=0 = 0' 

we can find that the Jacobian at the zeroth order solution has the 
form of: 

Jn 

h\ 



where 



^12 = - 



L^'Rf 



3R*/-A)L'- + W[A+R' 



hi 



L^R*^ 



{A + Rf) 



V [A + R%; 



-L^A 



(33a) 



(33b) 



The determinant of this Jacobian, given through: 
detj = -J12J21, 



(i) the "float" case, when the axis of symmetry of the rigid body 
is always perpendicular to the orbital plane, 

(ii) the "spoke" case, when the axis of symmetry lies along the 
relative position vector, 

(iii) the "arrow" case, when the axis of symmetry lies in the 
plane formed by tangential and normal to the orbital plane vectors. 

Now we may show that all these solutions are, in our terminology, 
locally central and Lagrangian, i.e., the point mass V, and the mass 
centre of the rigid body 2 move in one plane. 

Since for the description of the system there has been chosen 
the reference frame fixed at mass centre of the point body T, the 
potential of the body, besides the distance from the mass centre, 
depends also on the orientation of the rigid body. Following Ki- 
noshita, we denote a parameter, which describes the orientation, as 
V. In this designation, v is an angle under which the point body 7 
is "seen" from the mass centre of the rigid one S. Note an impor- 
tant fact that the "float" and "arrow" types of motion correspond 
to a case when t/v vanishes. Let us determine the physical meaning 
of that condition in terms of our parameterization. For that reason, 
we will use an expansion of the rigid body potential through se- 
ries in terms of the Legendre polynomials, the same as in |Kinoshita| 
( |I970[ |. This representation of the potential of 2 is the following: 



U ■ 



k=0 



(34) 



where v is the cosine of the angle between the symmetry axis of the 
body and the relative radius vector, and Pj; is the Legendre polyno- 
mial. For the "float" and the "arrow" types of motion, Uy vanishes. 
Hence, 



1 ^ 

f/v = - 1; A 

k=0 



■4(v)=0, 



for a certain value of v. 

In terms of variables used in our paper, the expansion of \34) 
reads as follows: 

Let us recall that a point is locally central if and only if 
pU3-R3Up=0. 
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Thus, we have: 

If f/v = then it means that the considered solution is locally cen- 
tral. Condition t/y = includes all possible locally central cases 
except of that ones when p vanishes. That is why the "arrow" and 
the "float" type relative equilibria found by IKinoshitaj j 1 970l l, un- 
der assumption that = 0, are locally central. Moreover, it is easy 
to see , the "spoke" solution is also locally central. In fact, for this 
solution we have p = t/p = 0. Thus, all solutions found by Kinoshita 
are locally central. 

The "float" equilibrium coincides exactly with our first so- 
lution as in both of these cases the axis of symmetry is always 
perpendicular to the orbital plane. The "arrow" and the "spoke" 
motions are just particular cases of our second solution. Indeed, the 
second solution, when SIt, vanishes, becomes the "spoke" motion, 
i.e., the axis of symmetry lies along the relative position vector; if 
we have = 0, then it coincides with the "arrow" motion. Quite 
surprisingly, in the paper of Kinoshita, we did not find any other 
cases corresponding to the inclined planar solution. It seems that 
investigating the "arrow" case, Kinoshita assumed for simplicity, 
that the body S is symmetric with respect to the equatorial plane. 
Moreover, then he used generally wrong implication saying that 
(/v = means v = 0. Let us notice that in our representation v = 
is equivalent to = 0. 

To summarize this Section, [Kinoshita] \ 1 970^ did not find the 
non-Lagrangian solutions because he implicitly assumed that one 
can always choose an inertial reference frame in such a way that 
the relative orbit lies in its (j;,)')-plane. In fact, this is equivalent to 
a priori assumption that all stationary solutions must be Lagrangian 
ones. 



7 CONCLUSIONS 

In this paper, we performed global, geometric analysis of the sta- 
tionary motions in the unrestricted problem of a point mass and 
an axially symmetric rigid body (the Kinoshita problem). Our aim 
was to determine all possible relative equilibria in this problem. 
We have shown that three types of stationary motions can be distin- 
guished. Two of them, the cylindrical precession and the inclined 
planar motion, are Lagrangian and are characterized by co-planar 
orbits of the mass centers of the bodies, whereas the third type of 
solutions, the conic precession, is non-Lagrangian. 

In the cylindrical precession the axis of symmetry of the rigid 
body is always perpendicular to the orbital plane. This type of sta- 
tionary motion was also found by | Kinoshita] jT970|, where it is 
called the "^oaf" motion. 

In the inclined planar precession, the axis of symmetry of the 
rigid body is inclined to the orbit. Two special cases of such kind of 
motion are found by |Kinoshita| ( [T97 0i . In the first case the axis of 
symmetry lies along the relative position vector, and in the second 
one lies in the plane formed by the tangent and the normal to the 
orbit. In Kinoshita ( 19701 these solutions are named as the "spoke" 
and the "arrow" motions, respectively. 

In the conic precession, the axis of symmetry lies in the plane 
formed by the relative position and the normal to the orbital plane 
vectors. Moreover, the point and the mass centre of the rigid body 
move in different parallel planes. This type of stationary motion is 
completely new in the problem analyzed in this work, although, as 
we noted above, such solutions have been constructed for specific 



non-symmetric potentials by |Abulnaga & Barkin| l |197^ . We found 
them thanks to a formulation of the problem through the minimal 
set of assumptions, basically implying only the form of the gravita- 
tional potential of the rigid body, besides the analysis of the prob- 
lem in the very basic settings of the Newtonian dynamics. 

We show that the determination of the stationary motions in 
the problem of Kinoshita may be reduced to solving at most two 
non-linear, algebraic equations. Our results can be applied to an 
arbitrary axially symmetric body, providing that an explicit form of 
its potential function is given. 

In our forthcoming paper l | Vereshchagin et al.|2009[ l, which is 
a direct continuation of this work, we perform the stability analysis 
of the relative equilibria found here, and we apply our approach to 
study the dynamics of a few specific models (i.e., choosing explicit 
form of the gravitational potential of the rigid body). 
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